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Abstract: This work addresses nonperturbative effects in botli matrix models and topologi- 
cal strings, and their relation with the large-order behavior of the expansion. We study 
instanton configurations in generic one-cut matrix models, obtaining explicit results for the 
one-instanton amplitude at both one and two loops. The holographic description of topological 
strings in terms of matrix models implies that our nonperturbative results also apply to topologi- 
cal strings on toric Calabi-Yau manifolds. This yields very precise predictions for the large-order 
behavior of the perturbative genus expansion, both in conventional matrix models and in topo- 
logical string theory. We test these predictions in detail in various examples, including the quartic 
matrix model, topological strings on the local curve, and Hurwitz theory. In all these cases we 
provide extensive numerical checks which heavily support our nonperturbative analytical results. 
Moreover, since all these models have a critical point describing two-dimensional gravity, we also 
obtain in this way the large-order asymptotics of the relevant solution to the Painleve I equation, 
including corrections in inverse genus. From a mathematical point of view, our results predict 
the large-genus asymptotics of simple Hurwitz numbers and of local Gromov-Witten invariants. 
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1. Introduction 



A well-known result in field theory states that the expansion of gauge theories has non- 
perturbative corrections which behave as e~^ ||l], 0. Physically, these corrections are due to 
instantons in the collective field theory which describes the large limit. In cases where the 
gauge theory has a string theory dual, they typically correspond to D-brane instanton effects. 
Based on rather general field theoretic arguments, one should expect that these e~^ effects are 
further related to the large-order behavior of the expansion, as is familiar in standard 
perturbation theory Q| . 

Perhaps the simplest class of large N gauge theories with string theory duals are matrix 
models. In spite of their apparent simplicity, matrix models hide a great deal of nontrivial 
information, as there are two different classes of string theories which can be described with 
these models. The first class of examples are the so-called noncritical or minimal string theories, 
defined as two-dimensional gravity coupled to conformal matter with central charge c < 1. To 
be precise, these theories are described by matrix models in the double-scaling limit, i.e., near 
critical points (see @ for an excellent review). The second class of examples which are described 
by matrix models are topological strings: as it was shown in [^, the genus expansion of the 
topological B-model, on certain noncompact Calabi-Yau (CY) backgrounds, is described by the 

expansion of a certain type of matrix models. Furthermore, it has been recently shown that 
the mirrors of toric manifolds can also be holographically described via matrix model technology 
1^ 1^. In this set of examples there is no need to go near a critical point in order to have a string 
dual; a generic CY background is here described by a matrix model off-criticality, whose couplings 
precisely correspond to the moduli of the CY in question. If one further tunes the matrix model 
couplings to a critical point, the matrix model will describe a particular CY background at, say, 
a given fixed value of the Kahler parameter. 

One of the main goals in our present work is to understand nonperturbative phenomena, 
as described by instantons, in both matrix models and topological strings (where, in the latter, 
we always have in mind their dual gauge theoretic description via matrix models). For the case 
of matrix models, the instanton configurations in the expansion have been identified long 
ago in terms of eigenvalue tunneling ||2|, ^, and they have been studied in great detail in 
the double-scaling limit. In [§, ^, David considered the action of an instanton configuration, 
which is obtained by analyzing the tunneling of a single matrix eigenvalue across the unstable 
effective potential — in this context this corresponds to a one-instanton effect; tunneling of several 
eigenvalues would correspond to multi-instanton effects. David explicitly showed that, near the 
critical point, this one-instanton action precisely agrees with the large-genus behavior of the 
free energy, which is in turn obtained from the matrix model in the double-scaling limit, via 
a solution to the so-called string equation. In the dual string theory, these effects were later 
identified as D-instanton effects [10, |l^] due to the so-called ZZ branes |12], and it was shown 
in ||To| , |l3| that a direct D-brane calculation reproduces the instanton action obtained from the 
double-scaled matrix model. This line of research thus made precise the connection between 
D-instantons in string theory and eigenvalue tunneling in the matrix model dual. 

Quantum fluctuations around this one-instanton configuration, again restricted to the double- 
scaling limit, were further analyzed in [^J], and more recently in |15, 16|, but the connection to 
the large-order behavior of perturbation theory was never explicitly addressed in any of those 
papers. In fact, it it surprising that to this date there has been no detailed study of instanton 
configurations in the matrix model per se, i.e., off-criticality, nor of their connection to the large- 
order behavior of the expansion. In [|l^, [l^ a general setting for this study was presented 
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but, unfortunately, the results of these papers, albeit written in terms of general matrix model 
data, are incorrect once we move away from the critical point. 

In this paper, we study in detail the perturbative expansion around a one-instanton configu- 
ration in a generic, one-cut matrix model. In particular, we shall give explicit formulae for both 
the one and two-loop contributions, and we shall write them in terms of geometric data which 
only depend on the spectral curve associated to the matrix model. This is a critical aspect of our 
analysis as it makes it possible to apply our results, not only in the realm of conventional one-cut 
matrix models, but also to more general theories which are defined by geometric constructions 
based on a spectral curve, as in [^]. In particular, and of special interest to us, this is the case 
of topological strings on certain toric backgrounds Q, and our general formulae make it possible 
to compute instanton effects in these models as well. 

Indeed, an important motivation for this paper is to use the dual matrix model description 
of topological strings on local CY manifolds as a nonperturbative definition, which then makes it 
possible to compute instanton effects in these theories for the first time. This is very similar to 
the nonperturbative holographic definition of noncritical string theories by double-scaled matrix 
models. By using this description, we deduce that the nonperturbative completion of the topolog- 
ical string theories considered in this paper includes an infinite number of nontrivial topological 
sectors, corresponding to the different instanton sectors of the matrix model. Geometrically, we 
interpret these nonperturbative effects as due to domain walls interpolating between D-brane 
configurations, as it had already been anticipated in 

A rather important aspect of all our nonperturbative computations is that they are testable 
via their connection to the large-order behavior of perturbation theory. Since we compute 
instanton effects up to (and including) two loops, we can determine the large-order behavior of the 
genus g free energy up to (and including) the 1/g correction. There are various examples where 
one can compute the expansion to high order, and by making use of standard numerical 
techniques which extract the asymptotic behavior of a perturbative series, we find an impressive 
agreement between the large-order numerical data and our theoretical instanton predictions. We 
will analyze in detail two types of examples. The first example concerns the standard hermitian 
quartic matrix model, studied for example in [18|. The second class of examples deals with 
topological string theory on local curves, which was extensively studied in ^. We shall 
confirm and improve the predictions of about the large-order behavior of these models, and 
we will also consider a special limit of topological strings on local curves which describes simple 
Hurwitz numbers (studied in [^0|). All of these models have a critical point, describing pure 2d 
gravity, which is controlled by the Painleve I equation. The double-scaling limit of our instanton 
calculations provides results for the large-order behavior of 2d gravity which refine those obtained 



1 21, ^ and agree with the analysis of the asymptotics in |23|. In fact, with the help of the Painleve 
I equation one can derive the full perturbative expansion around the one-instanton sector, and 
in this way we provide a further check of our explicit two-loop calculation. 

Mathematically, our results are highly nontrivial predictions for the asymptotics of the 1/A^ 
expansion of a one-cut matrix model, and they provide some clues concerning the analytic struc- 
ture of the total free energy of topological string theory, as a function of the string coupling 
constant. In the case of topological strings, our tests of large-order behavior provide a fur- 
ther check of the conjecture in as well as new conjectures about asymptotic properties of 
enumerative invariants that have not been explored so far. 

This paper is organized as follows. We begin in section 2 by presenting a short review of 
instanton effects and their connection to the large-order behavior of perturbation theory. We 
review a simple quantum mechanical example and further provide an extension of the main ideas 
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Figure 1: The potential for tlie quartic anharmonic oscillator. When g > the theory has an unstable 
vacuum at the origin, which decays via instanton tunneling. This vacuum gets stabilized when g < 0. 



to the expansion and string theory. These ideas are then explicitly applied in the analysis 
of one-instanton effects in matrix models in section 3. Here, we shall follow the general strategy 
put forward in but we shall both simplify and considerably improve their results. In 

particular, we shall present complete formulae for both the one-loop and the two-loop corrections 
around the one-instanton configuration, in generic one-cut matrix models. Applications of these 
results are then considered, starting in section 4 where we consider the quartic matrix model 
both off-criticality and in the double-scaling limit where it becomes pure 2d gravity. We further 
present numerical tests of the predictions given by the instanton calculation, by analyzing the 
large-order behavior of both the quartic matrix model and the Painleve I equation. We then 
proceed to consider applications in topological string theory. In section 5 we shall consider 
topological string theories on local curves, verifying and extending the predictions of Q, and 
we shall discuss the spacetime interpretation of the instanton effects in terms of domain walls. 
Then, in section 6, we analyze in detail the large-order behavior of the generating functionals for 
simple Hurwitz numbers as a further example of our formalism. In all cases, we find impressive 
agreement between theoretical and numerical results. A concluding section presents a list of open 
problems raised by our work. Finally, we also collect some explicit formulae for the free energies 
of both the quartic matrix model and Hurwitz theory, at high genera, in an appendix. 



2. Instantons and Large-Order Behavior 

In this section we shall review the connection between instantons and the large-order behavior 



of perturbation theory. Good references on this subject include |^, 24, 25, 26]. 
2.1 Field Theory Models 

Let us start by considering a quantum mechanical or field theoretical model which depends on 
a coupling constant, g, in such a way that for g > the theory has an unstable vacuum and 
that this vacuum gets stabilized for g < 0. A simple example of such a situation is the familiar 
quartic anharmonic oscillator with potential 

V = ^x^-gx\ (2.1) 

Due to the instability, there will be instanton solutions (sometimes called bounces in this context) 
which mediate the decay of the false vacuum. This is illustrated in Fig. As one analytically 
continues the coupling constant to the full complex plane, one finds that the partition function will 
have a branch cut along the real, positive g axis, with a discontinuity which is purely imaginary. 
In particular, one may write for the full partition function 

Z{g±ie) = Z^^\g)±^discZ{g), (2.2) 
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Figure 2: The complex plane for the functional integration. Here, and arc the rotated contours 
one needs to consider for 17 > 0. Their sum may be evaluated by the contribution of the saddle-point 
at the origin. Their difference is evaluated by the contribution of the sub-leading saddle-points, here 
denoted as Si and 62- 



defining both Z^^') and the discontinuity across the branch cut disc Z{g) = Z{g + ie) — Z{g — ie). 
A careful analysis of the physics of this problem, in the particular example of the anharmonic 
oscihator P, |2|, |2l, shows that is given by the path integral around the perturbative 



vacuum (or zero-instanton configuration), while the leading contribution to disc Z{g) turns out 
to be given by the path integral calculated around the one-instanton configuration, i.e. the 
instanton configuration with the lowest action in absolute value. We shall denote this path 
integral by Z^^\g). 

Let us be slightly more precise on this point. If we want the partition function to remain 
meaningful, as one performs the analytical continuation in the coupling constant from the stable 
to the unstable case, it is required that the contour of integration is also rotated, in a compensat- 
ing way [|^] {e.g., in the quartic oscillator as one continues —g to — 5exp(ibi7r) one must rotate x 
to X exp(=pi7r/4)). The rotated integration contours are illustrated in Fig. ^. What the analysis in 
|27, 28, ^ shows is that Z^^^ is computed as the integral over the sum of both contours, C"*" +C~ . 



In particular, if one is to compute the path integral in a saddle-point approximation, the contri- 
bution to arises from the saddle-point at the origin. On the other hand, the discontinuity 
disc Z{g) is computed on the difference of the two rotated contours, — . This immediately 
implies that the saddle-point at the origin cancels, between the two contours. One thus needs 
to consider the sub-leading saddle-points, which correspond to the one-instanton configuration. 
These sub-leading saddle-point contributions are also illustrated in Fig. |2|. In particular, notice 
that 

zW(g)~e-^/^ (2.3) 

and it is exponentially suppressed for small g as compared to Z^^^ . This is exactly as one should 
expect from the discussion above. If we now consider the free energy, defined by -F = log Z, we 
similarly have 

F{g ± ie) = fW(9) ± \disc F{g), (2.4) 
where fW(c/) = log Z^°\g) and 

disc F{g) = log = + . . . , (2.5) 
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at leading order in e '^/S" , "We will denote by 

- 1^ (...) 

the one-instanton contribution to the discontinuity. The zero-instanton sector has a perturbative 
expansion around g = given by 

+00 

fW(5) = J^afc/, (2.7) 

fc=0 

while the contribution from the one-instanton sector to the discontinuity disc-F(g) turns out to 
have an expansion of the form 



F«(5)=i<7-'e-^/^5^c„9". (2. 



n=0 

In this equation, A is the action of the single instanton, 6 is a characteristic exponent, and c„ is 
the (n + l)~loop contribution around the instanton configuration. 

If one now assumes analyticity of F{g) in the i^-plane, except for the branch cut along the 
positive real axis which we alluded to before, as well as some suitable conditions on the g ^ 00 
behavior, one can deduce the following relation between the coefficients of the perturbative 
expansion around the zero-instanton sector and the discontinuity across the cut 

Plugging the expansion for F^^^ (^]1) in the above formula (|2.9| ) we find an asymptotic expansion 
for large k, 

-, +00 

"fc ~ ^ E CnA-''-'+^Tik + b-n). (2.10) 

n=0 



This can be equivalently written as 



A 



Uk ~ — T{k + b) 

Ztt 



ciA C2A^ 
k + b-1 {k + b-2){k + b-l) 



(2.11) 



What one learns from this analysis is that the computation of the one-instanton partition func- 
tion, at one-loop, determines the leading order of the asymptotic expansion for the perturbative 
coefficients of the zero-instanton partition function, while higher-loop corrections yield the 1/k 
corrections. Notice that instanton configurations with an action A' > A (in particular, multi- 
instanton configurations with action nA, n > 2) give corrections to the asymptotics of Ofc which 
are exponentially suppressed in k, and will not be considered in here. The relation between a 
nonperturbative instanton computation and the large-order behavior of perturbation theory was 
first implemented by Bender and Wu in the case of the quartic anharmonic oscillator in quan- 
tum mechanics |2^]. They used the WKB method in order to perform a two-loop computation 
around the bounce, and thus obtain precise numerical values for cq and ci. Furthermore, they 
performed accurate numerical tests of their prediction (12.91) for the large-order behavior of the 



Ofc coefficients. Their results were later reproduced in path integral language |27] 



In this quantum mechanical example the analyticity conditions for the free energy can be 
justified rigorously (see |30| for a review). In more general situations (such as in quantum field 
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theory) one cannot justify these same assumptions; however the relation (|2.9| ) can be tested in a 
number of examples with surprising numerical precision (see, e.g., p4| , p^ ] for a review of these 
tests). 

2.2 The Expansion and String Theory 

The existence of a connection between instantons and large-order behavior has also been ad- 



dressed in the context of the expansion; for example in |31| where one considers vector 
models in low dimension. In the case of matrix models and their double-scaling limit, such a 
connection was used in |8|, ^ ^ in order to infer on the large-order behavior of pure 2d gravity, 
by computing the instanton action directly in the matrix model (see [Q] for a review). However, 
precise tests at one-loop or higher (the Cn coefficients in the expressions above) have not been 
performed to date, and we shall fill such a gap in the present work. In order to proceed to loop- 
level, one first needs a generalization of both the dispersion relation ( |2.9| ) and the expression for 
the perturbation theory asymptotics ( |2.10|) , to the present setting. 

We shall proceed in a heuristic way. Let us first consider the perturbative series in the 
zero-instanton sector of a closed string theory or its matrix model dual, 

+00 

F(^\gs) = J2F,{t)g'r'- (2-12) 
9=0 

In this equation the sum is over all genera, gs is the string coupling constant and t is the 't Hooft 
coupling t = gsN in the context of matrix models, or a geometric modulus in string theory. 
Observe that while in the previous case of the anharmonic oscillator one wanted to study the 
asymptotics of a standard numerical series, one now wants to address the asymptotics of a series of 



functions, naturally enlarging the complexity of the problem [31|. In order to have a perturbative 
series with standard structure, we shall consider instead 

T{gs) = glF{gs). (2.13) 

In this case, the one-instanton path integral yields a series of the form 

:F^'\z) = iz-^/2g-^ 2 (2.14) 

n=0 

where z = g1. This is an important feature distinguishing matrix models and string theory 
from field theory: the action of an instanton goes like and not as 1/z. Similarly, the 

perturbation series around the instanton sector is a series in powers of y^, and not a series in 
powers of z. As such, we may now write 

+00 

.FW(z) = ^F,(t)z^. (2.15) 

Our basic assumption is that a dispersion relation of the form ( ^^ ) holds in here, as it did in 
field theory. In this case, one finds 

= ^ £ - \ £ ^nA"^^-^""r(2, + 6 - n), (2.16) 

n=0 n=0 
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which may be exphcitly written as 



TT 



r(25 + 6)/ii 



1 + 



+ 



25 + 6-1 {2g + b-2){2g + b-l] 



+ ■ 



where we have introduced for later convenience 

^1 = Co, 



Ci 



= — , ^ > 1- 

Co 



(2.17) 



(2.18) 



The series inside the brackets in ( |2.17| ) must be understood as an asymptotic expansion in powers 
of 1/g, therefore up to two loops we can write it as 



IT 



T{2g + b) fii 



1 + ^ + 



(2.19) 



Justifying that the dispersion relation ( |2.9| ) holds in the present context is more delicate. The 
underlying reason is that g^ or l/N"^ appear naturally as coupling constants only in a collective 
field treatment of the problem (or, equivalently, in a formulation in terms of a closed string 
field theory). In spite of this, one could still present a heuristic derivation of ( |2.16| ) by making 
use of the Lipatov approach to the large-order behavior, and applying it within the context of 
collective/string field theory. In this approach one does not use the analyticity properties of the 
free energy, but instead performs a saddle-point evaluation in both field space and coupling space 
1 32]. Another heuristic derivation of ( |2.16| ) can be done by using Borel transforms Instead 
of trying to provide a more rigorous foundation for (|2.16| ) , we shall proceed to test it in various 
examples, also in the spirit of the many tests performed in field theory. 

In writing ( 2.17| ) we have implicitly assumed that there is a single instanton solution that 
contributes to the asymptotic behavior. In general there might be various instanton configura- 
tions in the system, with the same action in absolute value, and in this case F^^^ will denote the 
sum of all these contributions. For example, in the quartic matrix model, which we will analyze 
in section 4, due to the symmetry of the potential there are two instantons which contribute 
equally. It is also common to have complex instanton solutions which give complex conjugate 
contributions to F^^\ and in this case the asymptotic behavior of Fg is again obtained by adding 
their contributions [pSl . If we write 



A = \A\e'^^, 

the leading asymptotics will read in this case 



Ail 



|/ii|e^^w^ 



Fg- 



TT 



T{2g + h) 1^1 1 cos((25 + 6)0a + 9^ 



(2.20) 



(2.21) 



We shall also find examples of this situation in the models studied in this paper. 
2.3 Numerical Methods and Richardson Transforms 

The instanton computations we perform in this work yield predictions for the quantities A, b, 
/ii and iJ:2 appearing in ( p. 17 ) above. In order to test these predictions, one has to extract 
these quantities from the asymptotics of the sequence {Fg}g>o. However, computation of the 
amplitudes Fg is, in most cases, rather involved and therefore they will typically only be available 
at low genus, of order g < 20. This will also be the case for our examples, apart from 2d gravity 
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where the Painleve I equation allows for a computation to arbitrarily high genus. We shall 



therefore use a standard numerical technique known as Richardson extrapolation (see, e.g., [34|), 
in order to be able to extract the asymptotic behavior more accurately from the very first terms 
of the series. This method removes the first terms of the subleading tail and thus accelerates 
convergence towards the leading asymptotics. 

The basic idea of Richardson extrapolation is as follows. Given a sequence 

5(g)=so + - + ^ + ---, (2.22) 

9 9^ 

its Richardson transform is defined as 

A:>0 ^ ' 

This cancels the sub-leading terms in S{g) up to order g~^ . Indeed, one can show that if S{g) 
truncates at order g~^ , the Richardson transform gives exactly the leading term sq. 

The first quantity that one may extract from the sequence {Fg}g->Q, assuming it is of the form 
(2.17), is the instanton action. In order to apply the Richardson method, we need a sequence 



with large g asymptotics of the form ( ^.22 ). This is achieved by considering the sequence 



Fg+i 1 , 1 + 26 , 1 



Once A has been found, one can then simply extract the parameter h from the new sequence 
Finally, one obtains the coefficients /ii and ^2 from the sequences 



Ml 1 + ^ + O K (2-26) 



r(25 + 6) ^'V 25 W 

and 

whose asymptotics are already of the form ( |2.22| ), with leading terms and /U2) respectively. 
This is the basic picture behind most of our numerical work. 

The situation is slightly more complicated when we have to deal with two complex conjugate 
instantons. In this case, the ansatz for Fg is given by ( 2.21| ). If the absolute value of the instanton 



action is known, its phase 9a can be checked using the sequence 

I \A\^^-^Fg.,i2g + 6 - 2)(2, + 6 - 1) _ ^ ^^^^^^^^ 



{2g + h + l){2g + h)Fg Fg y y 

(2.28) 

3. Instanton Calculus in Matrix Models 

We shall now perform a more systematic implementation of the ideas discussed in the previous 
section, in the context of generic, one-cut matrix models. 
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3.1 Preliminary Results on Matrix Models 

We will consider a hermitian matrix model for an x matrix M, with generic potential V{M). 
We shall use the normalizations of |35|, so that the partition function will be defined by 

where the factor vol(C/(iV)) is the usual volume factor of the gauge group that arises after fixing 
the gauge. In terms of eigenvalues in the diagonal gauge, Zjy reads 



^^=m^j ndA.A^(A)e-i^^-^(^'\ (3.2) 



i=l 



where A(A) = ni<j(Aj ~ ^j) is the familiar Vandermonde determinant. The normalized free 
energy of the matrix model is then defined by 

F = log|g, (3.3) 
where is the partition function of the Gaussian matrix model, defined by the potential 



V{M) = \M'^. The free energy has a perturbative genus expansion 



F = Y,F,{t)gl^-\ (3.4) 

9=0 

where 

t = gsN (3.5) 

is the 't Hooft coupling. Another important set of quantities in a matrix model are the connected 
correlation functions 

H^/.(pi,...,pO = (Tr^— •••Tr^— \ , (3.6) 

where the subscript (c) means connected. These correlation functions are generating functions 
for multi-trace correlators of the form 

Wh{pu ...,Ph)=J2 n,+i ^ n,+i (Tr M"^ . . . Tr M"'^)(,) , (3.7) 
and they further have a gs expansion of the form 

WhiPl, . . . ,Ph) = ^gl'^^'-^Wg^hiPi, . . . ,Ph;t). (3.8) 

g=0 

In the literature one may find a great deal of work concerning the computation of the quantities 
Fg{t) and Wg^h{pi, ■ ■ ■ ,Ph'it), starting with the seminal work of and culminating in the recent 
formulation of |17]. In the following, we shall focus on the so-called one-cut matrix models and 
present some well known results which will be required at a later stage in our computation. 
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At large N the zero-instanton sector, or trivial saddle-point, of the matrix model is char- 
acterized by a density of eigenvalues p{X). In the one-cut case, this density has support on a 
single, connected interval C = [a,b] in the complex plane. This density is completely determined 
by the condition that the effective potential on an eigenvalue, 



V,s{X) = V{X)-2t J dAV(A')log|A-A'|, (3.9) 

has to be constant — at fixed 't Hooft coupling — on the interval C: 

V,s{\)=tm, AgC. (3.10) 
A quantity which is closely related to the density of eigenvalues is the resolvent, defined by 

iVo{p)= [ dX^. (3.11) 
J P- X 

Once the resolvent is known, the eigenvalue density follows as 

piX) = (^o(A + ie) -iOoiX- ie)) . (3.12) 
It turns out that the resolvent may be written as 

Mp) = Yti^'ip)-yip)), (3.13) 

where y{p) is a function which has a branch cut along C, called the spectral curve of the matrix 
model. It is explicitly given by 

y{p) = M{p)^{p-a){p-b), (3.14) 
where M{p), known as the moment function, is given by 

MM^/|im \^ . (3.15) 

Joo 27ri z-p ^{z -a){z- b) 

with the contour of integration being around the point at oo. The endpoints of the cut follow 
from the equations 

dz V'{z) 

2^i^-a)iz-b) 

^ (3.16) 

dz zV'(z) 

fc 2vri ^{z-a){z-b) 

There is also a useful formula for the moments of the function M{p), which are defined as 



= why. ^"^"^ 



k>l, (3.17) 

p=a,b 



given in terms of contour integrals 



^ 2Tn (^^_af+-2{z-b)-2 ' Jc^T^^ iz-a)Hz-b)^+-2 
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To make a long story short, it turns out that the quantities Fg{t) and Wg^hiPi, ■ ■ ■ ,Ph',t) can 
be computed in terms of the spectral curve alone. More precisely, knowledge of the endpoints of 
the cut, a and b, and of the moments ( 3.17 ), is all one needs in order to compute them. This 
was first made clear in [37| and later culminated in the geometric formalism of p8| , 3£, 17|. For 
example, one has for the genus-one free energy [^] 



F, = --log [M{a)M{b){a-b)^]. 



(3.19) 



The two and three-point correlators at genus zero are given by |40| 



Wo,2ip,q) 
Wo,3{p,q,r) 



1 



pq - k{P + + b) + ab 



^{p-q? \ V{p-a)ip-b)iq-a)iq-b) 
1 

8\/ {p -a){p- b){q - a){q - b){r - a){r - b) 



(3.20) 



+ 



1 



^M(a) {p - a){q - a){r - a) M(b) [p - b){q - h){r - b) 
while the one-point function at genus one is given by l3 



Wi,i{p) 



1 



2p + 5-3a M'{a) 



+ 



lQM{a){p -a)^{p-a){p-b) \{p - a){b - a) M{a) 
1 / 2p + a-36 M'{b)\ 



+ 



l6M{b){p - b)y/{p -a){p-b) \{p - h){a - b) M{b) J ' 
The only exceptions to this rule are the genus-zero free energy, -Fo(i)) which is given by 

Fo{t) = -\j d\p{\)V{\) - \t\{t\ (3.21) 
J c 



and the one-point function 



Woa{p) = tujoip). 



(3.22) 



It is clear from the expressions above that, for these two quantities, the spectral curve is not 
enough and one also needs to know the explicit form of the potential. As we shall soon unfold, 
the perturbative expansion around the one-instanton solution is again completely determined by 
the geometry of the spectral curve. 

In the one-instanton computation that we shall perform in the next section, we will also 
need some results about the derivatives with respect to t of various quantities that characterize 
the large solution. A result we need is (see Q) 



d{tLJo{p)) 



(3.23) 



dt ^(^p-a){p-b) 
together with the following derivatives, which follow from the defining relations ( 3.16| ) and ( 3.18 ), 

(3.24) 



da 
dt 



1 



a-b M{a) 



db 
di 



1 



b-a M{b) 
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Using these formulae one finds, 



dty{z) = , 

^J(z - a)(z - b) 

^ (3.25) 

BM(z)- ^ ( {z-b)M{z) {z-a)M{z) \ 

* ^ ' {z- a) {z-h)\{a- b)M{a) (b - a)M{b) ) ' 

as well as 

^ , 6 M'{a) 2 / M{a)\ 

{6.1b) 

Finally, we will also need derivatives of the free energies. One finds^ lH 

dtFo{t) = -tm = -V,Q{b), 
d^Foit) = -dtV,s{b) = 2 log (5 - a) - 2 log 4, 

while higher derivatives with respect to t follow from ( |3.24| ). 

3.2 The One Instanton Sector of the Matrix Model 

We shall now compute the "path" integral around the one-instanton configuration of the matrix 



(3.27) 



model. In this process, we will adopt the framework put forward in |15], but we shall use saddle- 
point technology rather than the approach based on orthogonal polynomials. Such a strategy has 
been considered before, as the approach of [O] was first rephrased in terms of the saddle-point 



perspective in [^] (for a third point of view, based on collective field theory, see [41]). At this 
stage, it is important to point out that our calculation will improve on the calculations in |15, 16] 
in three different ways. First of all, we shall fully exploit the saddle-point technology in order to 
present a much more succinct derivation of the final results. Secondly, we shall compute explicit 
formulae for the quantum expansion around the one-instanton solution up to two loops. Thirdly, 
and more importantly, we shall correct both the approach and the one-loop result in ]|l^, |l^ 
which, as they stand, are incorrect once one moves away from criticality. 

We thus consider a one-cut matrix model in which the effective potential has the form 
depicted in Fig. ^. It is constant along the cut C = [a,b], where there is a local, unstable 
minimum, and it has a local maximum at the point xq. The standard 1/A^ expansion is computed 
by considering the saddle-point configuration in which all of the N eigenvalues have support in 
the cut C. As was first pointed out in ^, ^, a /c-instanton configuration corresponds to a 
distinct saddle-point, in which N — k of the eigenvalues remain with support in the interval C, 
while k eigenvalues are placed at the local maximum xq (of course one still assumes that k <^ N). 

The matrix integral for the one-instanton sector (this is the case where a single eigenvalue 
sits at Xq) is given by p5|] 

N—l 

^-=mhwl dxe-i^^(^)/ ndA.A^(x,A„...,A._,)e-^S- (3.28) 



^Notice that 9tKiff(6) can be obtained by integrating dty{z) from — cxd to b, an integral which diverges logarith- 
mically. Sensible results are obtained [|l^ by always considering its regulated version, where one simply drops the 
divergent log z term. 
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M{xo) = 




Figure 3: The effective potential for the matrix model eigenvalues. 



where the first integral in x is over the nontrivial saddle-point contour, which we have denoted 
hy X G I, while the rest of the N — 1 eigenvalues are integrated around the standard saddle-point 
contour Iq. At this stage, it might be useful to have in mind Fig. |2| and the discussion concerning 
the quartic oscillator in section 2. The overall factor of N in front of the integral is a symmetry 
factor, counting the possible distinct ways of choosing one eigenvalue out of a set of N. One 
can easily write similar integrals for the /c-instanton contribution (see |jl5|) but we shall not 



m 



consider them in here. Taking all normalization factors into careful account, one finds for (3.28) 

4) = j;^^ {2^r-\N - 1)1 , £^dx (det(xl - M' 



x2\(0) - — V{x) 



(3.29) 

ix f{x). 

The notation in this equation is as follows. is the partition function evaluated around the 



standard saddle-point, and within the standard 1/N expansion. M' is an [N — 1) x [N — 1) 
hermitian matrix, and all of its eigenvalues are still integrated around the standard saddle-point. 
{0)f^ is the normalized vacuum expectation value of the gauge-invariant operator O, again 
computed around the standard saddle-point, 

lO)^^ - /AeXonf^idA.A^(A)0(A)e-^^-"-^^(^') 
/,,^„nf=idA.A2(A)e-^^-^(^') 



Finally, we have also defined 



fix) = (det(xl - Mrf^\ e-^^(^). (3.31) 

As we have seen in (^^), the one-instanton contribution to the free energy may be expressed in 
terms of the partition function, at leading order in e~^/^, by 

. ^(0) 

In the rest of this section we shall present a careful computation of this quantity. 
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In order to calculate the instanton contribution ( |3.32| ), we shall first compute /(x), as defined 
above in ( 3.31| ). Making use of the familiar relation det(2;l — M) = exp (tr ln(xl — M)) we obtain, 



(det(xl - Mf) = exp 



V — ((trln(xl-M))' 

Z — / el 



'(c) 



(3.33) 



which is written in terms of connected correlation functions (recall that the cumulant expansion 
precisely relates the generating functional of standard correlation functions to the generating 
functional of connected correlation functions as in this equality). The correlation functions 
appearing in ( p. 33 ) are nothing but integrated versions of the correlators in ( |3.6| ), evaluated 
at coincident points. Let us define 



Ag,h{x-t) 



-4ri(x; t) — ^ ^ 
fc=0 



dpi 



,Vh, 



Xl = --- = Xh=X 



m-2k+l 



(3.34) 



{n-2k + 1) 



n-2k+l 



(x; t), n> 1. 



In this notation, the general perturbative formula for the determinant follows as 

(det(xl - Mf) = exp ( ^ '-4„(x; t) ) , 

\n=0 / 

where An{x;t) is the n-loop contribution. We have, for example, 

Ao{x;t) = 2Ao,i(x;t), 
Ai{x;t) = 2Ao,2(.x;t), 
4 

A2ix;t) = -Ao^3{x;t) + 2Ai^i{x;t), 
2 

Asix; t) = -Ao4x; t) + 2^i,2(x; t). 



(3.35) 



(3.36) 



One observes that in order to compute the determinant at n-loops, one would require analytic 
expressions for the Wg^h with ((7, h) = (0, n + 1), (1, n — 1), (2, n — 3), . . . , (^, 1). Let us also point 
out that the integration constants involved in the integrations in ( |3.34| ) may be simply fixed by 
the large x expansion of the correlators. Indeed, we have the expansion 



((trln(xl-M))^)(,)= ^ 



(TV • • • Tr M"^ ) (,) ^-i 



(3.37) 



Next, we define the holomorphic effective potential, which combines the matrix model potential 
together with Ao{x; t), as 

Vh,es{x; t) = V{x) -2t j dpuJo{p) = V{x) - 2t j dp p{p) log(x - p), (3.38) 

which satisfies 

ViM^;t) = y{x) (3.39) 
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as well as 

ReVi,Mx;t) = Ve{iix), (3.40) 
where Ves{x) was earlier defined in ( |3.9D . Altogether, one finally has for the integrand 



(3.41) 



fix) = exp ( -ll4,efr(x;t') + V 5r'^n(3;; t') ) , 



where 



t' = gsiN-l)=t-gs. 



(3.42) 



This shift in the 't Hooft parameter is due to the fact that the correlation function involved in 
(3.31) is computed in a matrix model with — 1 eigenvalues (recall we removed one eigenvalue 
from the single-cut). Since we are computing the one-instanton contribution in the theory with 
N eigenvalues, we thus have to expand (3.41) around t. This gives further corrections in gs, 
which we make explicit as 



fix) = exp - g',-' ^ a,^H,efr(x; t) + E E 9s^''' Hi 

\ k=0 ' n=l k=0 ' / 



(3.43) 



We shall write this expression as 

fix) = exp 

where we define 

^x)^Y.9s~''^ni^ 



1 



Vh.six) + ^ix] 



(3.44) 



E 

n=l 



9s 



(_-\)n-l \k 

^ ^ d^v^M^)+Y,^^d'^Ar-kix) 



nl 



n=l 

One finds, for example, 

$i(x) = Aiix) +dtVh^esix), 



k=0 



(3.45) 



1 



C^2(X) = A2ix) - dtAlix) - - a,Vh,eff(x), 

$3(x) = Asix) - dtA2ix) + 1 d^Aiix) + ^ dfVi^Mx) 



(3.46) 



In expression ( 3.44| ) all quantities now depend on the standard 't Hooft parameter t for the model 
with N eigenvalues, and we have thus dropped the explicit dependence on t. The derivatives 
with respect to t can be performed by using the formulae we presented at the end of the last 
subsection. 

One may now proceed with the integration of fix), 



[ dx exp ( -—Vh^eaix) + ^ix)] 



(3.47) 



If we wish to evaluate this integral as a perturbative expansion around small string coupling, 
gs, we can do it using a saddle-point evaluation |15, The integration contour is over the 
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Figure 4: The spectral curve y{x) has a singular point at the nontrivial saddle a;o- 



nontrivial saddle characterizing the one-instanton sector, which is defined by the usual saddle- 
point requirement 

Vies{^o)=0 y{xo) = 0, (3.48) 

with xq located outside of the cut. If we use the explicit form of the spectral curve ( p.l4| ) we 
find the equivalent condition 

M(xo) = 0. (3.49) 

The saddle-point xq is typically a local maximum of the effective potential, as depicted in Fig. ^. 
Geometrically, the spectral curve is a curve of genus zero pinched at xq, as shown in Fig. ^. This 
was observed in |42| in the context of spectral curves for double-scaled matrix models, and their 



relation with minimal strings. Of course, it can happen that there is more than one solution to 
(|3l9| ). In this case, there will be various instantons and we will have to add up their contributions 



(the leading contribution arising from the instanton with the highest action, in absolute value). 
The calculation of ( 3.47| ) is now completely standard, and it reduces to Gaussian integrations. 



The result is 



l^^dxfix) = Jy^^^ exp (^-ll4,efr(xo) + $i(xo)) U + ^5"/nj , (3.50) 



where the /„ can be systematically computed in terms of the functions ^nix) and their deriva- 
tives, evaluated at the saddle-point xq, by making use of the Gaussian integral and the Gaussian 
moments. This is a long and tedious process, where one should be very careful with factors of gs- 
In particular, one splits the integrand into the standard Gaussian integrand plus the rest, where 
the rest should be power-series expanded in order to produce Gaussian moments. This process 
is source to some extra factors of Qs that must be properly considered. In any case, there are no 
conceptual difficulties in taking this calculation to arbitrary order. In order to find an explicit 
expression for the two-loop contribution to the one-instanton path integral, we shall later need 



2 1 7 5x^h,cfr(a;o) + <9^Vh,cfr(a::o)9a;$i(xo) \ H 



2(n%(^o)) ' ' J 24(k;^^(xo)) 



(3.51) 



Observe that the required evaluation of derivatives at xq, in the expression above, is a rather 
straightforward exercise as we are dealing in this case with rational functions. 
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The last ingredient needed to compute the one-instanton contribution is the quotient of 
partition functions in the expression for discF, ( p.32D . This quotient can be written in terms of 
the standard, perturbative free energies, since 



Z 



(0) 



z 



(0) 
N 



exp {F(i!) - F{t)) . 



(3.52) 



In the rest of this section both F{t) and Fg{t) shall denote the unnormalized free energies, i.e., 
F = log Z]\f. If one explicitly expands in gs, by both writing the above expression in terms of the 
standard 't Hooft parameter t alone, and further expanding the free energy in its perturbative 
genus expansion (|3.4D, it follows 



Z 



(0) 
N-l 



z 



(0) 
N 



L2J /_-.N„_2fc+l 



\n=Q 



k=0 



(n-2fc + 1)! 



Fk{t). 



(3.53) 



One has, for example. 



Go 
Gi 

G2 



-dtFoit), 

-^^d^^F,{t)-dtF,{t). 



(3.54) 



Putting together ( ^.50 ) and ( 3.53| ) above, we finally find that F^^^ has the structure 



1 Hi exp 




(3.55) 



n=l 



Collecting results above we obtain the following contributions to F^^\ up to two loops: 

A = Vh,efT(a:o) - Go{t), 



1 



V 2vrK:efr(^o) 

A*2 = /2 + G2it). 



exp ($i(xo) +ai(t)). 



(3.56) 



Let us now give explicit expressions for these quantities in terms of data associated to the 
spectral curve 

( CT) - First of ah, by using (|]5|), (|]2^ and ( |339| ) we find 



A = Vh,es{xo) -Vh,e{iib) = J dzy{z) 



(3.57) 



which is the instanton action (here, we use the fact that T4i,efr(&) = ^fr(^))- Notice that, as 
pointed out in this expression also has a geometric interpretation as the contour integral of 
the one-form y{z) dz, from the endpoint of the cut C to the singular point xq (recall Fig. Q). 

We next move to the one-loop contribution, and begin with the computation of ^i{x). One 
can find the result for ^0,2(2;; t) (which enters in the expression of ^1) simply by integrating the 
first formula in ( ^^Oj) 

- {a + b)/2 



^0,2(3;; t) = log 1 + 



y/{x - a){x - b) 



log 2. 



(3.58) 
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Using dm) one further finds, 

dtVh^esix) 



-4 log 



+ 4 log 2, 



and both these results together is all one requires to obtain 

<l>i(x) = — log (x — a)(x — b) 



(3.59) 



(3.60) 



Adding to $i(x) the result for Gi{t), which follows from (|3.27| ), it is simple to put all expressions 
together and obtain the contribution, /ii, of the one-loop fluctuations around the one-instanton 
configuration, 

(3.61) 



^ 27rM'(xo)[(xo-a)(xo-6) 



This formula is valid for any one-cut matrix model with an unstable potential. Notice that if xq 
is a local maximum of V^g (x), one will have that M'(xo) < 0, and hence ^ii will be real. It is also 
important to point out that our result (3.61) is different from the result obtained in [|l5|, |l6j. The 
reason is that, in these references, no distinction is made between correlation functions computed 
at t' and those computed at t. Correspondingly, the contribution of ( |3.5E ) is never taken into 
account. While this contribution vanishes at the critical point, it is non-zero for generic values 
of the parameters, making it crucial in order to obtain a generic result. In this paper we shall 
present substantial evidence that ( |3.6l| ) is the correct result, by using the connection to the 
large-order behavior of perturbation theory explained in the last section. 

The computation at two loops does not present any conceptual difficulty, but it is much more 
involved. One needs the explicit expressions 



1 



(a - by 



+ 



M(a)(x 
1 



.3 

a)2 



M{b){x-b)2 



12(a-6)2 \M{a) M{b) 



+ 



+ 



1 



(2(x - a) + (6 - a)) Vx-b {2{x - b) + {a - b)) ^/x 



(3.62) 



24(a - 6)2 
y/x — a — 



.3 

a) 2 



+ 



M{a) (x 
T f2M{a) + (a - b)M'{a) 



8(a - bf 



M2(a)V^ 



M{b) (x - 6)2 

2M(6) + (6 - a)M'{b) 
M2(6)^/^^ 



After very long but straightforward computations, one finally obtains the two-loop coefficient as 



1 



iJ'2 



(xo - b)M'{a) (xo - a)M'(6) 



M2( 



+ 



+ 



4(a - 6)y^(xo - a)(xo - 6) 

(xo - a)(xo - 6) f 8(xo - a) + 17(a - 6) 
12(^ 

5(M"(xo)) 



+ 



M2(6) 
8(xo - 6) + 17(6 



a) 



bf V {xo - afM{a) 
^ 3M'(xo)M(3)(xo) 35(2x0 



+ 



(xo - 6)2M(6) 
(a + 6)) M"(xo) 



24 (M'(xo))' ^{xo-a){xQ-b) 48 (M'(xo))^ ((xo - a)(xo - 6)) 
140 (2x0 - (a + 6))^ + 33(a - 6)2 



+ 



-+ 



(3.63) 



96M'(xo)((xo-a)(xo-6))^ 
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As one immediately realizes from the explicit expressions above, both fii and /i2 depend uniquely 
on data specified by the spectral curve. More precisely, they depend on the endpoints of the cut, a 
and b, the position of the saddle-point, xq, and on the moments of the function M(p), evaluated 
at a, b or xq. It is not hard to convince oneself that the rest of the coefficients fin in (3.55) 
must also share this property. This has two important consequences. First of all, it displays the 
universality of the results, in the sense that two matrix models which lead to the same spectral 
curve will also share the same discontinuity, disc F. In particular, since taking the double-scaling 
limit commutes with the geometric computation of the amplitudes, different models that lead 
to the same critical theory will also lead to the same one-instanton contribution at criticality 
1 15, 16 1 . Secondly, since the description of the B~model on mirrors of toric manifolds in ^ ^ only 
depends on the geometry of the spectral curve, we may also compute nonperturbative effects in 
these models by simple application of the formulae above for F^^^: one just has to apply them to 
the spectral curves described in |^, ^. Notice that in this paper we have restricted ourselves to 
the one-cut case and as such our formalism will only apply to the mirrors of local curves, worked 
out in B. 



4. Application I: Quartic Matrix Model and 2d Gravity 

Before proceeding towards the realm of topological string theory, we shall test our results in the 
case of a rather familiar matrix model, the quartic matrix model both off and at criticality. 

4.1 The Quartic Matrix Model 

The quartic matrix model is defined by the potential 

V{z) = ^z' + Xz\ (4.1) 

with A the quartic coupling constant. The properties of this model at large N were addressed long 
ago in |3€, 18]. The density of eigenvalues has support on the single cut C = [a,b] = [—2a, 2a], 
where a is a function of A and the 't Hooft parameter t, as 



1 



— ( -1 + Vl + 48At) . (4.2) 



24A 

The spectral curve follows as 

y{z) = M{z)^/z^ -4a'^, (4.3) 

with 

M{z) = l + 8Xa^ + 4Xz'^. (4.4) 
This function has two zeros which give two non-trivial saddle-points, namely itxo with 

xl = -^{l + 8Xa^). (4.5) 

These two saddle-points are evident in Fig. |5|, where we have displayed the effective potential for 
the quartic matrix model. If we wish to compare the large-order prediction of our formulae with 
the real behavior of the perturbation theory in this model, one is required to actually compute 
the free energies at high genera. The set-up for such a calculation was first described in |18| , 
but the calculation was only carried out in that paper up to genus g = 2. We have extended this 
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Figure 5: The effective potential Ves{x) for the quartic matrix model. There are two saddle-points 
located at xq and —xq. 



computation to genus g = W and we shall now review how to compute Fg in the quartic matrix 
model at large g. 

The calculation of the 1 /N expansion of the free energy in the quartic matrix model was set 
up in [|^] using the method of orthogonal polynomials. A review of such method and subsequent 
calculation would lead us too far apart from the main line of this work, so that in the following 
we restrict ourselves to presenting an algorithmic prescription to compute Fg which summarizes 



the results of The interested reader should consult the original reference |18| for full details. 



Also, for simplicity, we set t = 1 in the following and will follow the exact same conventions as 



in In particular, in this section our convention for the free energy, following [18|, is that 

F = —logZ. There are several components that make up the calculation of Fg. It starts with 
the so-called pre~string equation 

+ 4A {Rn + Rn^i + Rn+i)} = ng^, (4.6) 

for the coefficients Rn which determine the partition function in the orthogonal polynomial 
formalism. One then considers a continuous version of these coefficients, corresponding to a 
family of polynomials, r2s{x; A), which, in light of the pre-string equation ( [4. 61) , satisfy a simple 
algebraic, recursive relation. For s = 

ro(x;A) = ^ (-l + Vl + 48Ax) , (4.7) 
while for s > the pre-string equation yields the recursive expression 

r2„(x;A) + 2 J] ', ' \ = 0. (4.8) 

k+p=n ^ J 

In this way it is rather simple to compute the polynomials r2s(x;A) to very high s. These 
polynomials are crucial in other to find Fg. Indeed, the general formula for the total free energy 
is 11] 

glF{\) = - [ dx (1 - x) log H(x; A) + W(A)- 

p=i 



x=0 
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where the function A) is precisely built using the r2s{x; A) polynomials as 



.2s 



r2six; A) 



s=0 



In the expression above, B2p are Bernoulli numbers and TC^X) is the function 



(4.10) 



n{x) 



log 



dn e 



log 



e' 



(4.11) 



An expansion of ( [4. 91 ) in powers of gs then yields explicit expressions for Fg. Moreover, this is 
an algorithmic prescription of calculation, which may be simply implemented with a symbolic 
computation program. This calculation was carried out analytically up to 5 = 2 in p8[| and we 
have implemented it in a computer program, obtaining in this way explicit results up to 5 = 10. 
A partial list of our Fg can be found in the appendix. Here, let us just recall that |l^ conjectured 
that, for genus g > 2, the general structure should be of the form 



Fg{a^) 



'1 



a 



(2-a2 



with Vgia^) a polynomial in such that 



Vgia^ 



I) 



(45-3)! 



(4.12) 



2 • 629-1 gl{g - 1)!' 



(4.13) 



We have checked this conjecture up to genus g = 10 and further found that the polynomial 
Vgi^a^) is of order 85 — 4 in a^. 

4.2 Instanton Effects and Large Order Behavior 

Let us now present explicit formulae for the terms contributing to the one-instanton sector, in 
the quartic matrix model. As we did before, in the following we will set t = 1 for simplicity. The 
first thing to notice is that, since the potential is symmetric, there are two instanton solutions, 
corresponding to eigenvalue tunneling from C to the two saddles itxo (see Fig. ^) . Both instantons 
have the same action, which is computed via direct integration of the spectral curve 



A 



'4(l-a2 



2 log 



\/3\/-2 + q2 + V-2 - q2 +log4(l-a2), (4.14) 



and therefore contribute equally to the large-order behavior. The one-loop contribution ^ui can 
be easily obtained from the general formula we derived before, but it has an extra factor of 2 in 
order to account for the two instantons. It reads. 



1 



1 



a 



34 (2-a2)4 (2 + a2)4 



(4.15) 



After some tedious but straightforward analysis, one likewise obtains for the two-loop contribu- 
tion ^2 



1^2 



4^/3 



(2 - a2)2 (2 + Q 



2^1 



40 - I2a^ - 21a^ - lOa^ 



(4.16) 



- 22 - 



Figure 6: The sequence yl/Qg with Qg as defined in (2.24) and the corresponding Richardson transforms 
for the quartic matrix model, at fixed values A = —0.005 (left) and A = —0.01 (right). The prediction 
for the leading asymptotics is given by the instanton action A{X), shown as a straight line. The error for 
5 = 10 is 0.01% at A = -0.005, respectively 0.0047% at A = -0.01. 



Using these formulae, we see that 6 = —5/2 in (|2.14|) , and the asymptotics of Fg{\) is then given 
by 



(4.17) 



The goal is now to compare this "theoretical" large-order prediction with the actual, "exper- 
imental" behavior of the 1/A^ expansion, using the results we have obtained for the free energies 
Fg up to (7 = 10, in the quartic matrix model. We first focus on the range of values of A where 
the instanton action, as well as the Fg, are real. This is precisely the interval between A = and 
the critical point A = — ^ (we will come back to this critical point in the next subsection). In 
Fig. I^Fig. ^ we have displayed the asymptotic values of the instanton action as well as the one 
and two-loop results for the quartic potential. This is done at specific values of the coupling. 
The graphs include results extracted from the original sequence Fg (the uppermost sequence of 
data, colored in red), and its Richardson transforms (colored in orange, blue and green), along- 
side with the prediction from instanton calculus. In Fig. |^ we have plotted the asymptotic values 
of Hi and ^2-, obtained as a function of A from the third Richardson transform, divided by the 
corresponding prediction from instanton calculus. It is rather clear that this quotient is very 
close to 1, with a small error of roughly 0.1% over most of moduli space. The larger error found 
at A ~ is due to numerical difficulties related to the divergence of the instanton action in this 
region. Indeed, at very small A, the Richardson transformations converge too slowly to fall on 
a horizontal line at low genus — in this case, we would need higher-genus data to obtain better 
agreement with the predictions. In any case, the complete set of displayed numerical results 
strongly supports our analytical predictions. 

Let us now consider the range of moduli space where A > 0. In this region the amplitudes 
Fg are still real, as are the endpoints of the cut ib2a. However, the saddle-points xq given in 
( |4.5| ) now become purely imaginary and conjugate to each other. This implies that there are now 
four instanton solutions, corresponding to eigenvalues tunneling from both endpoints of the cut 
to both of the saddle-points, as depicted in Fig. The corresponding instanton actions are 
complex conjugate by a constant shift of ztivr. We therefore expect the leading asymptotics to 
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Figure 7: The sequence TTFgA'^^~^ /r(2g ~ |) and its Ricliardson transforms for tiie quartic matrix 
model, at fixed values A — —0.005 (left) and A = —0.01 (right). The prediction for the asymptotic value 
is the one-loop result /Lti (straight line). The error is 0.003% at A = —0.005, respectively 0.002% at 
A = -0.01. 



-0.20 



-0.25 



-0.5 



-0.6 




-0.7 - y 



Figure 8: The sequence ( ^.27| ) and its Richardson transforms for the quartic matrix model, at A = —0.005 
(left) and A = —0.01 (right). The prediction for the leading asymptotics is given by the two-loop result 
fi2. The error is 0.05% at A = -0.005, respectively 0.016% at A = -0.01. 



be of the form ( 2.21| ), implying that 



^Fg\A\ 



(4.18) 



where 6a and 9^-^ have been defined in ( 2.2C| ). This is indeed the case, as one can see from 
Fig. ^ showing the quotient in the left hand side of ( |4.18 ) together with the prediction for 



2 cos [{2g — |)^A + at two positive values A = 0.004 and A = 3. 
4.3 2d Gravity and the Painleve I Equation 

A rather well-known result (see for an excellent review) is that the quartic matrix model has 
a critical point at 

(4.19) 



A. 



1 

'48' 



At this critical value of A, the saddles itxo collide with the two endpoints of the cut ib2a. One 
may further use the matrix model near this point in order to define two-dimensional gravity by 
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Figure 9: The left figure shows the asymptotic value of irFgA^^^'i /r{2g — |) for the quartic matrix 
model, as extracted as a function of A by the third Richardson transform, divided by the analytic prediction 
fii. The figure on the right shows the analogous quotient for /Z2. For A < —0.004 the error is always less 
than 0.06%. 




Figure 10: This figure shows the instanton effects for the quartic matrix model at positive coupling A. 
The endpoints of the cut, a and b, are real while the two saddle-points of the effective potential are purely 
imaginary and complex conjugate to each other. There are two pairs of complex conjugate instantons, 
corresponding to eigenvalues tunneling from either end of the cut to the saddles xq and Xq. 



means of a double-scaling limit. In this specific limit, one takes 

A ^ Ac, gs ^ 0, (4.20) 

in such a way that the variable 

z = -^iX-X,)g;^/' (4.21) 

is kept fixed. In this limit it follows that the total, perturbative free energy of the matrix model 
becomes the free energy of pure 2d gravity 

F{gs,X)^FA,{z). (4.22) 

Furthermore, in this limit, the pre-string equation of the quartic matrix model (j^^) precisely 
becomes the Painleve I equation 

- -u" = z, (4.23) 
3 

governing the specific heat of the model 

u{z) = -F:[^{z). (4.24) 
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Figure 11: The sequence 7rFg|Aps~i/ (|^i|r(2g— |)) for the quartic matrix model, together with the 
prediction for the leading asymptotics 2 cos ((2g — |)6'^(A) + 9^^ (A)) (thin black Hne), at A = 0.004 (left), 
respectively A = 3 (right). At the highest depicted values of g the error is of the order of 2% (A — 0.004), 
respectively 5% (A = 3). 

These results may be used to obtain the perturbative expansion of Ffi^{z), at any given order. It 
turns out that the free energy obtained in this way is actually doubled, since it gets contributions 
from the two collisions at ^xq. This is of course due to the symmetry of the potential, which we 
have discussed before. In order to remove the doubling it is enough to change the normalization 
of the quantities appearing above, by 

z 2tz, u 2^u, Fds 2Fds. (4.25) 
Proceeding in this way one is led to the Painleve I equation with the normalization 

V? - ^u" = z, (4.26) 



while the double-scaled free energy still satisfies ( 4.24| ). The perturbative expansion of the specific 
heat has the form 



+ 0O 

u{z)=z'^Y,UgZ~^3/2^ (427) 

9=0 

SO that the Painleve I equation becomes equivalent to the following difference equation for the 
coefficients Ug 

25(9-1)2-1 1^ 
^3 = 4^ Ug-i - -2^utug-t, uo = l. (4.28) 

£=1 

The coefficients a^, which appear in the perturbative expansion of the double-scaled free energy 
as 

F,s{z) = -y/" - ^log . + Y^agZ~'^^-')'\ (4.29) 
can then be obtained from Ug through the simple relation 
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As a result one finds, for example, 



^ . N 4 5 1 7 

-ftisl-^) = logZ H 

^ ^ 15 48 5760 



Z 2 + 



245 
331776 



z~^ + 



(4.31) 



We are now in a position where we may obtain a prediction for the asymptotics of the coefficients 
of this series, a^, by simply evaluating the expressions we obtained for the quartic matrix model 
near the critical point, and taking into account the change of normalization in (4.25). In this 
way we find 

«^ ' (4.32) 



5 

Z4 . 



Moreover, for /ii we obtain 



while at two loops we get the result 



9s 



/9s 1^1 



8 • 34^ 



5 

Z 8, 



9s 



37 
64\/3 



_ 5 
Z 4. 



(4.33) 



(4.34) 



Altogether, this means that the one-instanton contribution to the double-scaled free energy, up 
to two-loop order, is 



F. 



(1) 



ds 



8 • 34^ 



z 8 exp 




37 
64^/3 



Z 4 + 



(4.35) 



The result for the one-loop coefficient, fii, was first obtained by David in |14] and later re-derived 
in |15]. We have obtained /i2 directly from an instanton computation in the matrix model, but 
we may also verify our result by computing the one-loop instanton expansion directly from the 
Painleve I equation. This expansion has been studied in detail in [23|, where it has been used 
to analyze the asymptotics of the perturbative answer. The calculation of this expansion goes 
as follows. As noticed in |22, ^, the discontinuity of the double-scaled free energy ( 4.35| ) can be 
computed by linearizing the string equation ( 4.26| ) around the perturbative, asymptotic solution. 
If we denote 

(4.36) 



€{Z 



discu{z), F£) = -e"(z), 



one finds the linear and homogeneous differential equation 

e"{z) - 12uo{z)e{z) = 0, 

where 



uo{z) = 

It is easy to solve ( 4.37| ) at 



1 

Z2 



1 



1 



-Z 2 



49 



1225 _i5 

:Z 2 + 



48 4608 55296 

oo after "peeling off' the exponential piece. 



(4.37) 



(4.38) 



e{z) = cz 8 exp 



8\/3 5 

Z4 



(4.39) 
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The overall coefficient c cannot be deduced from the differential equation ( [4.37 ) due to its homo- 
geneity, but the ek can be easily found in terms of the coefficients of the asymptotic expansion 
of uq. One finds, for the very first terms. 



eiz) 



cz s exp 




64^/3 



5 

"4 + 



75 
8192' 



341329 
23592960^/3' 



15 



. (4.40) 



Of course, the coefficient c may still be fixed with the explicit result for the one-loop coefficient 
/ii . Assembling all together, one finds the full perturbative expansion of the free energy around 
the one-instanton configuration. 



F. 



(1) 



1 



ds 



3 

34 



z 8 exp 



TT 




37 

64^/3 
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'2 + 



6433 
24576' 



12741169 
23592960^3' 



15 



(4.41) 

We can now use ( 4.35| ), together with ( |2.17| ), in order to obtain a prediction concerning the 
large-order behavior of the perturbative coefficients of the double-scaled free energy, as 



16^/30 / 25 Y 
125vri V192j 



37 



3927 



3618769 



80g 12800^2 15360000c/': 



+ 



(4.42) 



Since one can compute these coefficients up to very large order, by using the Painleve I equation. 



we can now perform truly precise tests of some of our proposals. In Fig. 12 we show numerical 



checks for both the one and two-loop predictions, up to genus 400. Indeed both leading and 
subleading asymptotics of the coefficients ag clearly agree, to a very high degree of precision, 
with our prediction ( 4.42| ). This also leads us to an important point. It is sometimes stated in 
the literature, e.g., |15, |T^, that the one-instanton amplitude (4.35) cannot be deduced from the 



Painleve I equation. The reason for this assertion is simply that the linearized equation ( 4.37 ) 
for e does not allow the calculation of ^i. But it is clear, in view of the connection between 
large-order behavior and instanton effects, that there is a more subtle relation between the one- 
instanton amplitude and the perturbative result. In fact, one could instead have derived this 
amplitude from the asymptotics of the coefficients ag, themselves derived from Painleve I. It is 
easy to see that from the difference equation ( [4.28| ) one may obtain 



192 j 



(4.43) 



a result which at leading order precisely agrees with ( [4.421) . A careful study of the difference 



equation (l4.28| ) beyond ( 4.43| ) conffims indeed the result ( 4.42| ) for the asymptotics of Ug, and 



in particular makes possible to extract the one-instanton amplitude directly from large order. 



5. Application II: Topological Strings on Local Curves 

We shall now proceed into the realm of topological string theory, beginning with the case of 
topological strings on local curves. 

5.1 Topological Strings and Matrix Models 

Topological strings are defined as twisted N = 2 sigma-models coupled to 2d gravity, and 
they provide a vast generalization of noncritical/minimal strings (where one considers conformal 
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Figure 12: The left figure shows the sequence agA^^^^^^^^^ /r(2(7 — 5/2) for 2d gravity alongside with its 
Richardson transforms, up to g = 400, clearly converging to the one-loop prediction 1/(8 • 3'^/'*7r'^/^) — 
0.009847. The right figure shows the modified sequence ( ^.27 ) and its Richardson transforms for 
converging towards the two-loop result — 37/(64\/3) 
genus is of order 10~^%. 



9' 



-0.33378, again up to g = 400. The error at this 



matter with c < 1 coupled to 2d gravity). As it is well known, topological strings come in two 
types, called the A-model and the B-model, which are related by mirror symmetry. The A-model 
provides a physical formulation of Gromov-Witten theory, while the B-model is deeply related 
to the theory of deformation of complex structures. In the last years it has become apparent 
that topological strings on CY threefolds share many of the remarkable properties of noncritical 
string theories, like integrability, but their geometrical structures are much richer. 

It is natural to ask if topological strings on CY threefolds have a matrix model description, 
at least in some cases. In Dijkgraaf and Vafa showed that the B-model, on certain CY 
geometries, is indeed equivalent to a matrix model (see |35] for a review of this development). 
This class of target geometries is of the form 

uv = H{X,Y), (5.1) 

where H{X, Y) is a polynomial in X,Y G C, and u,v £ C The nontrivial information about 
this geometry turns out to be encoded in the Riemann surface T, described by H{X, Y) = 0. An 
important insight of the analysis of Dijkgraaf and Vafa is that the spectral curve of the matrix 
model is precisely the Riemann surface S, providing in this way a beautiful example in which 
the master field of the 1/A^ expansion generates the target geometry of a string theory. 

More recently, the correspondence between matrix models and topological strings was ex- 
tended to toric CY threefolds 0]. This class of examples is very interesting since (in contrast 
to the geometries considered in Q) they have mirror geometries. These geometries are CY three- 
folds described by an equation of the form ( |5.1| ), but where the variables X,Y now belong to 
C*. The proposal of ^ is to regard these geometries as spectral curves of a matrix model. 
The open and closed string amplitudes of the B-model are then computed by the 1/A^ expansion 
associated to the spectral curve. Notice that one does not need to specify the matrix integral 



in order to compute these amplitudes; using the results of |37, 17] it is enough to specify the 
spectral curve in order to compute the 1/A^ expansion. 

In the previous section we have computed one-instanton effects in one-cut matrix models 
in terms of data associated to the spectral curve. We can then use the correspondence of |6|, ^ 
to apply our results to topological string theories described by this class of matrix models. The 
restriction to the one-cut case still leaves a rather general class of CY backgrounds to explore, 
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the so-called local curves. A special limit of the theory of local curves gives the theory of simple 
Hurwitz numbers studied for example in [20|, which will be addressed in the next section. 



5.2 Topological Strings on Local Curves 

Local curves are toric CY manifolds of the form 

Xp = 0{p - 2) © 0{-p) ^ P\ p E Z. (5.2) 

Topological string theory on Xp has received a lot of recent attention (see, e.g., [^] and references 
therein). As explained in the A-model on Xp has to be defined equivariantly, and the most 
natural choice (the equivariant CY case) corresponds to the antidiagonal action on the bundle 
(we refer the reader to |44| for further details). Of more interest to us in the present work is that 
the free energies at genus g on this geometry, Fg ^{t), depend on a single complexified Kahler 
parameter t, associated to the complexified area ofP^ They can be computed in both the A 
and B-models. 

In the A-model, the total partition function is given by 

+ 00 

Zx, = exp {Fx, {gs ,t)) , Fx, {gs , t) = ^ gs'~^F^^ (t) , (5.3) 

9=0 

and near t ^ oo, F^^ {t) has the expansion 



F^^{t) = Y,Ng,a^-''\ (5.4) 

d=l 



where Ng^d are the Gromov-Witten invariants of the CY manifold Xp at genus g and degree d. 

The total partition function Zx, can be computed as a sum over partitions, by making use 
of the topological vertex formalism as described in |4^. In order to write the explicit resulting 
formula, we first have to introduce some notation. To begin with, define the q-number [n] as 

[n] = g'^/^ _ g = e^'^ (5.5) 

A representation, R, of U{oo) is encoded by a Young tableau, labeled by the lengths of its rows 
{li}. The quantity 

i{R) = ^h (5.6) 

i 

is the total number of boxes in the tableau. Another important quantity associated to a given 
tableau is 

KR = ^li{li-2i + l). (5.7) 

i 

We finally introduce the quantity 

Wr = q-'^-'' n TTT^T^' (^•8) 

with hook(n) the hook-length. With all this notation at hand, we may finally write the explicit 
expression for the topological string partition function on Xp, which is given by 

Zx, = Y.^R^R'^^'"^^^'''^Q'^''^^ Q = (-l)Pe-*, (5.9) 

R 
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where ii* denotes the transposed Young tableau {i.e., the tableau where we have exchanged the 
rows with the columns). 

Although ( |5.9| ) gives an all-genus expression, it is effectively an expansion in powers of Q. 

— X 

In order to obtain an expression for each Fg ''(t) to all orders in Q, one usually appeals to mirror 

symmetry and the B-model. However, standard techniques of mirror symmetry do not work well 
when applied to local curves. The Riemann surface encoding the mirror geometry for local curves 
was proposed in Q based on the direct analysis of the sum over partitions presented in |19], and 
later on some aspects of this mirror construction where confirmed from a more mathematical 



point of view |46]. The B-model geometry is encoded in the spectral curve 



2/(A) 



tanh 



-1 



V(A - a)(A - b) 



ptanh 



-1 



V(A-a)(A-b) 
A + Vab 



(5.10) 



which has genus zero. Although this curve is not algebraic, it is easy to see that when written 
in terms of the variables X = X,Y = one obtains an algebraic equation for the C* variables 



X,Y, which leads to the mirror CY threefold (|5.l| ). The advantage of writing the curve in the 
nonalgebraic form (|5.10| ) is that, as explained in ^, one can apply verbatim the standard 
matrix model technology that we use in this paper. The curve ( p. 10 ) may also be written in 
the form ( ^.14| ), with a moment function M(A) which has various nontrivial zeroes where the 
spectral curve is singular. The endpoints of the cut, a and 6, are given by 



(5.11) 



where C is related to Q by the mirror map |19, ^] 

Q = (l-C)-P(P-2)^. 



(5.12) 



Xr, 



It was further conjectured in M that the free energies Fg ^ (t) can be obtained as the standard 



genus g free energies of a matrix model with spectral curve ( |5.10 ). And it was conjectured in 
|19] that, for g > 2, these free energies may be written as 



where 



Cc 



5(9-1) 



1=1 



1 



(5.13) 



(5.14) 



{p-iy 

is a critical point of the model. In fact, at this point, a zero xq of M{\) collides with the endpoint 



of the cut 6, and we are left with a critical theory in the universality class of pure 2d gravity jl^ . 
If one further takes the double-scaling limit. 



9s 



0, z fixed, 



where 



4(1-Cc)3' 

then the total free energy ( |5.3D becomes the free energy of pure 2d gravity. 



(5.15) 
(5.16) 
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5.3 Instanton Effects and Large Order Behavior 

In M the matrix model description, based on the spectral curve (5.10), was used to study 



nonperturbative effects in this topological string theory. The spectral curve ( 5.10 ) has a nontrivial 
saddle xq, which is the solution to 

M(xo) = 0. 



For the cases p 
by 



and 



3 and p = 4 the relevant solutions have been determined in H 



Aab 



Xo 



Xo 



2^h 



p = 4. 



(5.17) 
they are given 

(5.18) 
(5.19) 



In it was argued that this saddle controls the large-order behavior of F^''{t), at any value of 
t. We shall now show that this is indeed the case, and that the one and two-loop results pi,2 
computed in terms of the spectral curve ( |5.1[l| ) control the subleading large g asymptotics. The 
instanton action for an eigenvalue tunneling from b to xq has already been computed in it is 
given by the rather formidable expression 



where 

F{x) 



2Li 



log(/i(x))(^log(/i(:E))-21og(l + 
2/i(x) 



A{Q) = F{xo)-Fia), 



2fi{x) 



(5.20) 



2Li2 



2fi{x) 



+ log 1 + 



2/i(x) 
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log log x- 
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2/2(x) 
\{^ + Vb)'^J ■ 2 



) -log(] 



+ Vbf 



(5.21) 



^^^^^^ + 2pLi2 + ?(logx)2 +plogiV^+Vbflog X, 



and 



fi{x) = V(ar^a)(ar^ + 



a + b 



(5.22) 



f2{x) = y/{x - a){x - b) + X + Vab. 



In these expressions a and b are the endpoints of the cut as usual, given in ( [5.11 ). 
The one and two-loop coefficients are again given by the general expressions 

i^M) and 



(3.63), as in the previous section. We shall now compare the analytic results to the large-order 
behavior of the perturbation series for the case of the local curve (p = 3). The leading and 
subleading asymptotic behavior of Fg should be given by the same structure found for the quartic 
matrix model ( [4.17| ). As explained in section we can independently test the predictions 
for the instanton action, as well as the one and the two-loop results, by applying Richardson 
transformations to the modified sequences (2.24)-( [2.27] ). Notice that all these quantities depend 
on the B-model modulus C. For simplicity, we shall restrict our analysis to the range < C < 
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Figure 13: The sequence y/l/Qg and the corresponding Richardson transforms for the local curve X3, 
at fixed values ( = 0.24 (left) and ^ = 0.15 (right). The leading asymptotics are predicted to be given by 
the instanton action A{(), shown as a straight line. The error for the available degree g — 8 is 0.014% at 
C = 0.24, respectively 0.025% at C = 0.15. 



Cc = |) where the endpoints of the cut, as well as the instanton action, are real. Fig. ^ shows 
the inverse square root of the sequence Qg in ( ^.24 ) and its first three Richardson transforms. 



at two specific values of The straight line is the prediction for the instanton action, A. As is 
evident from the plot, and even though we only use data up to genus g = 8, the third Richardson 
transform already falls on the straight line. The mismatch between numerical extrapolation 
and the prediction is of order 0.02%. Analogously, we may check the one and two-loop results. 
In Fig. 14 and Fig. 15 we plot the modified sequences, ( 2.26 ) and ( |2.27| ), together with the 



corresponding Richardson transforms, again at two fixed values of the Kahler modulus. As 



explained in section 2.3, the predictions for their leading asymptotics are /ii, respectively /i2, 
which are shown in the figures as straight lines. Again, this is confirmed by the Richardson 
transforms, clearly converging to the prediction from instanton calculus. The error in here is of 
order 1%. 



Similar graphs can be produced at any other point in moduli space. Fig. 16 shows the 
asymptotic value of the instanton action, as approximated by the third Richardson transform, 
divided by the corresponding analytical prediction, and plotted as a function of the modulus 
over < ^ < 1/4. This quotient is indeed very close to one, as it should be from our discussion. 
Similarly, in Fig. |l^ we plot the asymptotic results for fii and //2, divided by the corresponding 
analytic predictions, as functions of C- Notice that while the agreement is excellent over most of 
moduli space, as one approaches C ~ the deviation from the predicted value increases. This 
is again due to the divergence of the instanton action at this particular point of moduli space. 
Indeed in this region, the Richardson transforms converge too slowly to fall on one line, at low 
genus g < 10. In order to obtain full agreement one would need higher-genus data, which is out 
of our scope in this paper. 

We have performed similar checks of our predictions for the local curve X4, also obtaining 
agreement to very high precision, and further strengthening our analytical results. 

5.4 Spacetime Interpretation of the Instanton Effects 



As we have seen in ( 3.57] ), the instanton action can be computed as a contour integral from 



the endpoint of the cut to the saddle xq. This contour integral measures the potential difference 
between the cut C and xq. When the spectral curve corresponds to a double-scaled matrix model, 
this instanton action should correspond to the disk amplitude for a D-instanton in noncritical 
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Figure 14: The sequence nFgA^^^^ /r(2.g — |) and its Richardson transforms for the focal curve X3, at 
fixed values C = 0.24 (left) and C, — 0.15 (right). The prediction for the asymptotic value is the one-loop 
result fii, shown as a straight line. The error is 0.49% at C = 0.24, respectively 0.58% at C = 0.15. 
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Figure 16: The asymptotic value of y^l/Qg for the local curve X3 as extracted from the third Richardson 
transform as a function of divided by the analytic prediction for the instanton action. For C > 0.05, 
the error is always less than 0.03%. 



string theory. These D-instanton configurations have been identified in terms of ZZ branes, and it 
has also been verified that indeed the matrix model computation agrees with the disk amplitude 
for a ZZ brane [|l3|. Equivalently, the ZZ disk amplitude can be calculated as the difference 
between the disk amplitudes for two FZZT branes located, respectively, at the branch cut of the 
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Figure 17: The left figure shows /ii for the local curve as extracted from the perturbative series using 
the third Richardson transform of the sequence ( 2.26D , divided by the corresponding analytical prediction, 
and plotted over the range < C < 1/4. Similarly, the second figure shows the asymptotic result for /i2 as 
obtained from the perturbative series using (2.27), again divided by the corresponding analytic prediction. 
The typical error is about 1.5%. 



curve and at the pinched point of the curve. It turns out that, for topological string theory on 
local curves, there is a similar interpretation of the instanton action in terms of D-branes, as 
well as a spacetime interpretation in terms of domain walls. 

The natural branes for the A-model on a toric CY manifold are the Harvey-Lawson branes, 
first studied in this context in |47|. The mirrors of these branes are just points in the spectral 
curve of the B-model. Two branes located at points zo and zi define an interpolating domain 
wall in the underlying type II theory. The tension of this domain wall is given by the difference 
of D-brane superpotentials |47| 

W{zi)-W{zo)= rdzyiz). (5.23) 

Jzo 

When Zo and zi correspond, respectively, to the endpoint of the cut and the saddle xq, ( 5.23 ) 
is exactly the instanton action computed in ( 3.571 ). The connection between instanton actions 
in the matrix model and tensions of domain walls was already made in |^] for the backgrounds 
considered therein. At the same time, (Hill) can be regarded as the difference between two disk 
amplitudes for D-branes located at zi and zq. We then see that the role of FZZT branes in 
noncritical string theory is played by the Harvey-Lawson branes in topological string theory on 
local CY threefolds. Indeed, it can be easily seen ^ that, in the case of local curves, the toric 
branes become FZZT branes near the critical point describing 2d gravity. On the other hand, the 
saddle xq that we have been considering (and which leads to an extremum of the superpotential) 
gives a topological string analogue of the ZZ brane. 

A more invariant way of writing ( 5.23| ), by taking into account the full six-dimensional 
geometry of the CY, is 

A = [ n, r=[Ci-Co], (5.24) 



where il. is the holomorphic (3, 0) form on the CY, and F is a three-cycle interpolating between 
the two-cycles Co,i associated to zq^i in the full geometry. This is indeed the general form for 
disk amplitudes of B-branes presented in |^] . 

It is interesting to notice that usually the nonperturbative effects due to B-branes considered 
in the literature involve the hypermultiplet moduli, since a B-brane supported on a curve will 
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couple to the Kahler form, and not to Q ||4^, 50 1. This type of D-instanton effects (which in 



some cases can be computed exactly cannot however be related to the large-order behavior 
of the topological string amplitudes, which depend on the vector multiplet moduli. On the other 
hand, domain walls interpolating between two B-branes can couple to 17 and therefore have the 
right structure to control the large-order behavior of topological string perturbation theory. In 
this paper we have checked this for a restricted class of toric geometries, but we expect this fact 
to be true in the more general case, for an appropriate choice of the domain wall. 

6. Application III: Hurwitz Theory 

We finally proceed to our last example, Hurwitz theory. 
6.1 Hurwitz Theory 

Hurwitz theory studies branched covers of Riemann surfaces. Here, we shall restrict ourselves 
to the coverings of a sphere P"^ (the "target") by surfaces of genus g (the "worldsheets" ) . The 
covering maps will be restricted to have only simple branch points. The number of disconnected 
coverings of degree d with these topological characteristics is counted by the so-called simple 
Hurwitz number, which we denote by ^^^^(1*^). It can be computed, in classical Hurwitz theory, 
in terms of representation theory of the symmetric group: 

2 



e{R)=d 

Here the sum is over Young tableaux R, with a fixed number of boxes i{R) equal to the degree 
d, and dn is the dimension of R regarded as a representation of the symmetric group 5"^. The 



quantity kr was defined in (|5.7D . 

We can now define the total partition function of Hurwitz theory as a generating functional 
for simple Hurwitz numbers, 

where Q = e~*^ and gn can be regarded as formal parameters keeping track of the degree and 
the genus, respectively. This partition function can be written as 

The free energy log describes connected, simple Hurwitz numbers Hg^^{l'^)', 



g>0 d>0 

and it has the genus expansion 



(6.4) 



F%H,tH) = Y.9H~^^fiQH)- (6.5) 

3=0 



-36- 



This theory is in fact a topological string theory in disguise. It can be realized as a special limit 
of the type-A theory on local curves Xp with Kahler parameter t that we studied in the previous 
section [^], namely the limit 



OO, t ^ OO, Qs ^ 0, 

while the new parameters gn and tn, which are defined by 

gH=pgs, e-*« = (-i)Ve-*, 



(6.6) 



(6.7) 



are kept fixed. As in the case of the theory on local curves, there is a B-model mirror to this 
theory. Its natural coordinate x is related to the A~model coordinate Q = e~*^ by the mirror 
map 



Q, 



(6.8) 



which can indeed be understood as an appropriate limit of ( 5.12| ) for p oo [19|. The inverse 
mirror map is provided by Lambert's W function |52], 



X 



-W{-Q) = J2 



k=l 



k\ 



(6.9) 



which has convergence radius Qc = e ^ or x = 1- The large-radius region corresponds to Q — > 

(6.10) 



(and also to x ~^ 0). The spectral curve characterizing the B-model is of the form 
y{h) = 2tanh~^ 



,Vi<^-h){b-h) 



2h-{a + b) 
where the endpoints of the cut are given by 



V(a-/i)(6-/i), 



i\2 



1\2 



(i-x^) 



(6.11) 



The above spectral curve can also be read from the saddle-point description of the sum over 
partitions ( |6.3| ) given in |53, 20 1. 

Hurwitz theory has been extensively studied in the mathematical literature, and these studies 
have unveiled interesting properties. As shown in |54[ |, the higher-genus free energies F^{Q), 
when expressed in terms of the mirror coordinate x, have a very simple structure, namely 



6 



sx' 



+ x, 



F?{x)=-^{log{l-x)+X 



(6.12) 



F^ix) 



Pgix) 



33-3 
i=2 



(1 - x)5(9-l) 

Moreover, the polynomials Pg{x) have the property 

Pg{l)=A3-^ag, g>2 



9>2. 



(6.13) 



where ag is the genus g free energy of 2d gravity appearing in ( 4.29] ). Therefore, in the double- 
scaling limit 

X^l, gn^O, g-\l-xf = AKl, (6.14) 
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the total free energy of Hurwitz theory becomes (|4.2S| ) 

F''{gH,tH)^FM, (6.15) 

and one recovers 2d gravity at the critical point. This was first pointed out at genus zero in ||2^ 
and then established at all genera in [ 19 1 , using the results of 1 54 1 . 

Another interesting result concerning Hurwitz theory was obtained in [55|, where the total 
free energy was shown to satisfy the Toda equation, 

exp {F^{gH,tH + gH) + F^{gH,tH - gu) - 2F^{gHM)) = g%e'd^^F^{gH ,tH)- (6.16) 



This equation is the analogue for this model of the pre-string equation (|4.6| ) for the quartic 
matrix model. One can directly derive from ( |6.16 ) that the double-scaled specific heat satisfies 
the Painleve I equation ( 4.26| ), providing in this way yet another derivation of the result in |19|. 
We have used the Toda equation to compute Hurwitz amplitudes up to genus 16, and some of 
these results are presented in appendix A. 2. 

6.2 Instanton EfTects and Large Order Behavior 

Let us now turn to the computation of the one-instanton quantities. From the curve ( 6.10| ) we 
find the moment function. 



M{h) 



: tanh 



-1 



^ia-h){b-h) 
where the nontrivial saddle-point is defined by 

M{ho) 



, Via-h)ib-h) 
2h-{a + b) 



or 



^/{a- ho){b- ho) 



tanh 



0, 



V(q- ho){b- hp) 
2ho - (a + 6) 



This equation can be written in a simpler way by defining w as 



ho 



w 



4y^cosh^(-) + (1- Vx)^. 



In terms of these variables, equation ( |6.19 ) simply reads 



w 



sinh(t(;) 



(6.17) 

(6.18) 
(6.19) 

(6.20) 

(6.21) 



Even though we cannot solve analytically for /io(x)) we can solve ( 6.19| ) to find /io(x) near x = 0, 1 
as a power series. Near the critical point x = 1) it is easy to see that /io(x) h^-s a Taylor series 
expansions in powers of = 1 — x 



(6.22) 



Near x = Oi power series solution is more complicated. At leading order it is easy to find 
that ^ 

w log(x) + log(- log(x)), (6.23) 
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which yields 



^o(x) ~ - log x + 21og(-log x), 



(6.24) 



The corrections to the leading asymptotics ( 6.23| ) can be obtained following a method exposed, 
for example, in |5£]. The full solution can be written as 



where f is a power series 



in the variables 



log(x)' 



w = log(- log(— )) - - log(x) + V, 



j,k,Tn 



log(-log(x)) 
T = ^ — — , ^ 



X 



log(x) Vlog(x:) 

The coefficients Cjkm can be explicitly written as 

dz e(J-i)^z'=+^(-l)'" {j + k + m)\ 



(6.25) 
(6.26) 

(6.27) 



Cjkm 



+ 



The instanton action 



27ri (e-^ - l)i+k+m+i j]k\m\ 

dz eP''z''{-l)"' {j + k + m-iy. 
27ri (e-^ - l)i+'=+'^ j\{k - l)!m! 

Aix) = / dhyih) 



+ 



(6.28) 



can now be computed explicitly as a function of ho as 

(a - 6)2 



Mx) = {b-a) (7 cosh ^(7) - VT^^ 



7 



VT^ — 1 — cosh ^(7) 



where 



7 



b — a 



i2hoix)-a-b). 



(6.29) 



(6.30) 



(6.31) 



Using the above results for the behavior of ho near x = 0, 1, we can also find the behavior of the 
instanton action near these points. At the critical point, one finds 



9H 5 



(6.32) 



where k is the double-scaled variable introduced in ( |6.14 ). Of course, this is the expected 
universal, double-scaled result of (4.32). Near x = 0, we find 



^(x) ~ ^(^og x)^ X ^ 0. 



(6.33) 



We may now compare our predictions with the numerical asymptotics of Fg. The moduli 
space of X can be divided into the six segments shown in Fig. 18, and we have tested our 
predictions in each of them. They have the following characteristics: 
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Figure 18: The moduli space of Hurwitz theory. 



The simplest case to study is the real interval < x < 1- Here, there is one single instanton 
with real action, corresponding to an eigenvalue tunneling from b to the saddle xq on the 
right of the cut, and all Fg are also real. 

As X moves to the right of [0, 1], beyond the critical point at x = !> the only solutions to 
M{h) = are located inside the cut and the instanton action becomes purely imaginary 
while /ii_2 remain real. The Fg oscillate in sign. 

It turns out that there is no systematic difference between the regions I- VI away from the 
real axis. The instanton action as well as ^1^2 and of course Fg are generically complex, in 
spite of which our predictions continue to hold. 

When X lies on the negative real line, the endpoints of the cut move away from the real 
axis and become complex conjugate. There are now two saddle-point solutions, xq and Xq, 
complex conjugate to each other, and accordingly two instanton solutions with conjugate 
actions, one corresponding to an eigenvalue tunneling from b to xq and another from a = b* 
to Xq, as shown in Fig. Therefore the Fg are real, with asymptotics of the form ( 2.21 ) 



involving a cosine. Notice that this is very similar to a mechanism for the local curve, first 
observed in M]. 



As before, the one and two-loop coefficients are given by ( |3.61|) and ( 3.63 ), evaluated for the 



moment function (6.17). The saddle-point solution has to be evaluated numerically. The instan- 



ton action, as well as /xi and //2) are well-defined over the whole complex plane of the modulus 
X- Fig. |2^ and Fig. 21 show the inverse square root of the sequence Qg in ( |2.24| ) and of the 



corresponding Richardson transforms, alongside with the prediction of instanton calculus for the 
instanton action, at values of the modulus x = 0.5, x = 1-5 + i, and x = —1 — 0.5i. In Fig. ^ we 
compare the sequence TrFg\A\'^^~^^'^ /{T{2g — |)|/ii|) for Hurwitz theory, together with the predic- 
tion 2 cos ((2(7 — ^) Oa + (^fii), at X = —0.5 and x = —3. Fig. ^ and Fig. 24 show the modified 



sequences, ( p. 261 ) and ( p. 27] ), with leading asymptotics given by the one and two-loop fluctua- 
tions around the one-instanton configuration, together with the analytic prediction. Indeed, the 
agreement is again quite spectacular. 

7. Conclusions and Outlook 

In this paper we have extended classical results on the connection between instanton effects and 
the large-order behavior of perturbation theory to the realm of general, one-cut matrix models 
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Figure 20: The left figure shows the sequence (Re{Qg ^))^^^ for Hurwitz theory, together with its Richard- 
son transforms. The straight line shows the corresponding prediction (Re(A^))^/^, at x = 0.5. On the 
right, the same for x = 1.5 -I- i. The available degree is 5 = 16, the error is 4 x 10~^% at x = 0.5, and 
7 X 10-6% at X = 1.5 -t-i. 



and topological strings. After having performed a detailed one-instanton computation, up to two 
loops, we have tested this connection in both the standard quartic matrix model off-criticality and 
in its double-scaled limit, 2d gravity. Combining our results with the holographic matrix model 
descriptions of topological string theory on toric backgrounds , we have further provided a 
computation of nonperturbative effects in certain topological string models, and we have verified 
in detail that they precisely capture the large-order behavior of the perturbative amplitudes. It 
is important to point out that the precise agreement we find up to two loops strongly supports 
our claim that we have identified important nonperturbative effects in these models. This is also 
a strong check of the instanton/large-order correspondence proposed in this paper, as well as of 
the proposal of |^, |^ for describing topological string theories on toric backgrounds in terms of 
dual matrix models. 

From the mathematical point of view, we have presented precise conjectures for the large- 
order behavior of the Fg amplitudes, which are well-defined objects in all the cases we have 
considered. The simplest case is of course the proposed asymptotics ( 4.42[) , for the free energy 
of pure gravity, but the asymptotics of Hurwitz theory (i.e., the Hurwitz numbers) should also 
be of mathematical interest as they may provide new insight into this enumerative problem. 

Our work raises various problems and, at the same time, suggests various venues for future 
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Figure 21: On the left, the sequences (Re{Qg with its Richardson transforms and the prediction 

(Re(A^))^/^, at X = —1 — 0.5i (straight line). On the right, we show the same for the imaginary parts. 
The errors at g = 16 are 0.01% and 0.08%, respectively. 
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Figure 22: The sequence TrFg\A\'^3 ^/'^ /{T{2g — |)|Aii|) for Hurwitz theory, together with the prediction 



-3 (right). The error at genus 16 is of order 3%. 
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Figure 23: The left figure shows ■nA^^~iFg/T{2g — |) for Hurwitz theory, and its Richardson transforms, 
at X = 0.5. The leading asymptotics are predicted by /xi, shown as a straight line. On the right, we plot 
the analogous sequence ( 2.27 ), together with the expected leading asymptotic value /i2 (straight line). 
The error at 5 = 16 is 0.009% for m, and 0.012% for ^2. 



research, to which we hope to return in the near future. Let us thus finish this paper by listing 
some of them. 



One model we have not studied in here is topological string theory on the resolved conifold. 



1.000 



1.010 r 



Figure 24: The left figure shows /ii as a function of x for Hurwitz t heor y, as extracted from the pertur- 
bative series using the third Richardson transform of the sequence (2.26), divided by the corresponding 
analytical prediction. Similarly, the second figure shows the asymptotic result for ^2 as obtained from the 
perturbative series using (2.27), again divided by the corresponding analytic prediction. For x > 0.01, 
the error is of order 0.1%. 



It is known (see, for example, |35] and references therein) that this theory can be described 
by a Hermitian matrix model with a potential of the form (log x)^, which has a global 
minimum at x = 1 and no saddles. In principle, the way to address the large-order 
behavior in such cases is to deform the model, in order to obtain an unstable potential with 
a calculable one-instanton amplitude. The original potential is then recovered by analytic 
continuation [33, B^]. However, we have not found a suitable deformation which makes it 



possible to find the large-order behavior for this potential. In fact, we would face the exact 
same problems if we were to address the large-order behavior of the perturbation series for 
the ground-state energy of a (log x)^ potential in quantum mechanics. This is a situation 
which, to the best of our knowledge, has also not been addressed in the literature. 

Our work may be straightforwardly generalized to more complicated matrix models. For 
example, one could consider two-matrix models. Instanton effects in two-matrix models 



have been computed at leading order in |57], and near the critical point in |55], but it 



would be interesting to have exact results beyond leading order and off-criticality. Also, 
a more geometric formulation of the instanton contribution computed in this paper, along 
the lines of the approach in [0], would be desirable. Although our expressions only depend 
on the form of the spectral curve, they are only suitable in principle for a genus-zero curve 
written in the form ( 3.14| ). 



In this work we have restricted ourselves to the analysis of the one-instanton sector, but 
the nonperturbative completion of the theory involves A;-instanton sectors corresponding 
to the tunneling of k eigenvalues. One can easily extend the framework of ||T^, adopted 
in here, to compute these effects. Results for the two-instanton sector near the critical 
point have been obtained in [^9|. It would be interesting to make a detailed analysis of the 
multi-instanton sectors off-criticality. 

A possibly interesting check of our one and two-loop results for 2d gravity involves re- 
deriving them directly in the continuum Liouville theory, where instanton effects are de- 
scribed by D-brane instantons with ZZ boundary conditions [10, 11, 12, 13]. 
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It would also be interesting to rigorously verify some of the proposals for the asymptotic 
behaviors that we have put forward, beyond the numerical tests already performed in this 
paper. It seems very likely that the asymptotics of pure gravity can be established via a 
detailed analysis of the difference equation ( [4.28 ). One may also envisage that the two-loop 



instanton computation could be checked using the methods developed in [SO, 62| 



More generally, it would be rather interesting to provide a more rigorous foundation for 
the validity of the dispersion relation ( |2.9D , in both cases of matrix models and string 
theory. We should add however that, to the best of our knowledge, this kind of relation 
has only been properly justified in quantum mechanical models. In most examples of field 
theory such a relation is simply assumed to be true, and later tested a posteriori by explicit 
computations [^]. In spite of this, there might be some hope to further develop this line 
of research in the specific case of matrix models. 

Nonperturbative effects of order e~^ can be also found in models defined by sums over 
partitions, such as the case of two-dimensional Yang-Mills l63, These effects have 



been used in holographic descriptions of topological string theory It would thus 

be interesting to see if there is any relation between this description and the one we have 
proposed (and carefully tested) in terms of matrix models. 

• The topological string models we have studied in this paper are not very conventional, 
since they correspond to toric diagrams with intersecting lines, and this is refiected in the 
fact that their spectral curve is pinched. One could smooth out these models by resolving 
the singularity, obtaining in this way a spectral curve of genus-one with two cuts. In 
the context of noncritical string theory, this process is interpreted as adding ZZ branes to 



the background [42|. It would be very interesting to see if this leads to some geometric 



transition from the local curve backgrounds to other topological string backgrounds. 

For us, the most pressing problem is to extend the analysis in this paper to more complicated 
topological string models. The first case to address is that of matrix models with multiple 
cuts. In fact, the general multi-cut model with fixed filling fractions can be regarded as 
a matrix model in a generic, fixed multi-instanton sector, and by studying nearby filling 
fractions one obtains a general framework to address multi-instanton effects in matrix 
models. Some aspects of this framework were discussed in [^^, albeit in a different context. 
Once the multi-cut case is understood in detail one could use the philosophy of this paper 
to extend the results to topological string theory. This would allow for computation of 
nonperturbative effects on new and interesting toric backgrounds, such as the case of local 
P^, and would also be an important step in further strengthening our understanding of 
nonperturbative effects in topological string theory. It might even give precious hints for 
the future study of these effects in compact backgrounds. 

Although we have used the matrix model description of |^, |^] to compute nonperturbative 
effects in topological string theory, we have not provided a full nonperturbative definition 
of these models in this work. The reason is that the effects we have described in here 
only depend on the geometry of the spectral curve, and thus it was not necessary to write 
any topological string theory partition function explicitly as a matrix integral. In fact, it 
might happen that there is more than one way to do this, since two matrix models with 
different potentials and different finite N partition functions might nevertheless have the 
same spectral curve, and therefore the same expansion and one-instanton amplitudes. 
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As such, it would be important to go beyond the calculation of instanton effects and 
provide a full nonperturbative, holographic definition of topological string theories on toric 
backgrounds, in terms of matrix integrals. 
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A. Explicit Higher— Genus Formulae 

In this appendix we present some explicit expressions for free energies at high genera, in both 
the quartic matrix model and Hurwitz theory. This is just a partial list of our results, as most 
formulae quickly become too intricate to put in print. In spite of this, we hope these explicit 
expressions may be of future interest (and, as far as we know, have never been computed before). 

A.l Quartic Matrix Model 

As we have reviewed in section 4, an algorithm for computing free energies in the quartic matrix 
model was put forward in |18], and we have applied it up to genus g = 10. Here, we present a 
partial list of our final results. In [|^], the quartic free energies were computed up to genus two, 
with the result (here t = 1) 

Foia') = -\ log (a^) - ^ (l - a') (9 - a^) , (A.l) 
Fi(a2) = -Llog(2-a2), (A.2) 



v3 



F,{a^) = 1 ji-^ (82 + 21a^ - 3a^) . (A.3) 

It was further conjectured that, for genus g > 2, the general structure of the free energies should 
be of the form 

with Vg{a^) a polynomial in such that 

•pr^2 ^^ 1 (4g-3)! 

=') = 2:6^^!(^- ^^-^^ 



Using the exact same procedure as in |]18[, we have extended the analysis up to genus ten, 
verifying both conjectures above. In particular, we have obtained at genus three 

V2,{a^) = --^^ (17260 + o? (-32704 + Oa^ (-325 + OSa^ - 15a^ + a^))) , (A.6) 

which can be explicitly compared to another genus three calculation performed in [^, |6^, with 
both results in complete agreement. At genus four, we obtained 

V4{a^) = — (-1421392 + (12438536 + (-13719796 + 270^ (-15694+ 

38880 

+5810a2 - 14560"^ + 238a^ - 23a^ + a^°)))) , (A.7) 



-45 - 



and at genus five 



1 



85536 



-383964880 + (-1573981616 + (7592114712+ 



+q2 (-6114807776 + Sla^ (-781725 + 3268110^ - 1019610"^ + 23535a^+ 
-3915Q* + 445a^° - 31a^2 + a^^))))) ■ (A. 

Finally, at genus six the free energy follows from the polynomial 
1 



n(a") 



139728961867968 + -369974786833952+ 
79606800 ^ ^ 

(-955888270184512 + Sa^ (1037832523698416 + (-662581722466844+ 



+55971a^ (-39761282 + 17910398a^ - 63711120"^ + 1787698a'' 
+65901a^° - 8214a^2 ^ ^Qa^^ - 39a^^ + a^^)))))) ■ 



392007a* + 



(A.9) 



Although we have extended this calculation up to genus ten, the expressions quickly get too 
messy and little illuminating, and as such we shall not display any further polynomials. Our 
results further allow us to conclude that the polynomial Vg{oP') is of order 3(7 — 4 in a^. One final 
consistency check concerns the case of a^ = 2, corresponding to the critical point of the quartic 
model. In this situation it must be the case that 



p^(a2 = 2) = (-l)3 25(^'-i): 



"5' 



(A.IO) 



where Ug are the coefficients appearing in the expansion of the double-scaled free energy (obtained 
from Painleve I in the 1/3 normalization; see section 4 for details). Again, our results pass the 
test. 

A. 2 Hurwitz Theory 

If we expand ( 6.16| ) in qh, and use that [^, [l^ 



X 



E 

k=l 



k\ 



-tuk 



(A.ll) 



we obtain the recursion relation 



\l>l k>2 l>0 ' ^ 



2g 
9h 



(A.12) 

where we keep the coefficient of in the right hand side. If we combine this recursion with the 
general form of Hurwitz numbers ( 6.12| ), we can obtain explicit expressions for the polynomials 
appearing in ( |6.12 ) up to high genus. The first few are. 



2KX) 240 1440' 



X 



6 



+ 



53x' 



1741x^ ^ 137x^ 



+ 



X 



1008 10080 362880 181440 80640 



(A.13) 
(A.14) 
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P( )- ^\ '^Q^^^' 42419x^ 87739xQ 280603x^ 109x' 1291x' 
^^^^ 1440 604800 1209600 2177280 17418240 53760 21772800 7257600' ^ ' ^ 

^ 17387x^^ 67289x^° 44696593x'^ 193701347^^ 37315313^^ 

^'•^^ ~ 1056 ^ 665280 ^ 345600 ^ 79833600 ^ 273715200 ^ 91238400 ~^ 
8679559x^ 2295119^' 1525901x' 23x' 

82114560 205286400 3832012800 7603200 958003200' ^ ' ^ 
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